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We analyze the entanglement properties of spins (qubits) close to the boundary of spin chains in 
the vicinity of a quantum critical point and show that the concurrence at the boundary is significantly 
different from the one of bulk spins. We also discuss the von Neumann entropy of dissipative envi- 
ronments in the vicinity of a (boundary) critical point, such as two Ising-coupled Kondo-impurities 
or the dissipative two-level system. Our results indicate that the entanglement (concurrence and/or 
von Neumann entropy) changes abruptly at the point where coherent quantum oscillations cease 
to exist. The phase transition modifies significantly less the entanglement if no symmetry breaking 
, field is applied and we argue that this might be a general property of the entanglement of dissipa- 

^ ' five systems. We finally analyze the entanglement of an harmonic chain between the two ends as 

function of the system size. 
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I. INTRODUCTION 



, Coherence, decoherence and the measurement process are longstanding problems of quantum mechanics since 
' they mark the fundamental difference to classical systems. They have gained increasing importance in the context 
^ ] of quantum computing because the operation of a quantum computer requires a careful control of the interaction 

■ between the system and its macroscopic environment. The resulting entanglement between the system's degrees of 
• ' freedom and the reservoir has been a recurrent topic since the formulation of quantum mechanics, as it is relevant to 

. the analysis of the measurement process [U, [3, Q • 

' Theoretical research on Macroscopic Quantum Tunneling lead, among other results, to the (re)formulation of a 
I , canonical model for the analysis of a quantum system interacting with a macroscopic environment, the so called 

■ Caldeira-Leggett model [3], initially introduced by Feynman and Vernonfs'l. It can be shown that this canonical 
Ch model describes correctly the low energy features of a system which, in the classical limit, undergoes Ohmic dissipation 

(linear friction). It can be extended to systems with more complicated, non- linear, dissipative properties, usually called 
sub-Ohmic and super-Ohmic, see below [1, 0I- 

In relation to the ongoing research on entanglement, a recent interesting development is the analysis of the con- 
currence of spin-chains like the transverse Ising model and the XY model. It was found that the derivative of the 
^1 , concurrence obeys universal scaling relations close to the quantum critical point and eventually diverges at the tran- 

■ sition [1, Also other models which exhibit a quantum phase transition were subsequently investigated in this 
' direction, as e.g. the Lipkin-Meshkov-Glick model [T]|. 

y—{ ' Originally, the concurrence as measure of entanglement was introduced by Wooters due to its accessibility. 
I Alternatively, the von Neumann entropy of macroscopic (contiguous) subsystems can be used [l^. A non-local 
«j . measure of entanglement was employed in the study of the Aflleck-Kennedy-Lieb-Tasaki (AKLT) model [l^. [Tsj. 
^>«^ ' In this paper, we will first discuss the effect of boundaries of the Ising model on the entanglement properties using 
, the concurrence as measure of entanglement. We will then discuss a model which exhibits a boundary phase transition, 
L| ' i.e., two Ising spins which are coupled to two Kondo impurities. This model can be mapped onto the spin-boson model 
. ^ ' and the concurrence can be computed which was originally defined for the two Ising spins at the boundary [l^ . We 
' will then discuss various dissipative systems and c omp ute the von Neumann entropy, focusing the discussion on the 
Jh , cross-over from coherent to incoherent oscillations |29| . 
. . . . The von Neumann entropy is a more general information measure than the concurrence since the latter can only be 
defined for two spin-1/2 systems. The former can further be generalized to a measure which relates non-contiguous 
sub-systems which shall be done in the third part of this paper in the context of an harmonic chain. We note by 
passing that the concurrence is an essentially local measure which yields zero for all spin pairs which arc not nearest 
or next-nearest neighbors. 

We close the introduction with some general remarks. The models studied here, i.e., also the spin chains, can be 
interpreted as quantum systems characterized by a small number of degrees of freedom coupled to a macroscopic 
reservoir. These models show a crossover between different regimes, or even exhibit a quantum critical point. As 
this behavior is induced by the presence of a reservoir with a large number of degrees of freedom, they can also be 
considered as a model of dephasing and loss of quantum coherence. It is worth noting that there is a close connection 
between models describing impurities coupled to a reservoir, and strongly correlated systems near a quantum critical 
point, as evidenced by Dynamical Mean Field Theory In the limit of large coordination, the properties of 
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an homogeneous system can be reduced to those of an impurity interacting with an appropriately chosen reservoir. 
Hence, in the hmit of large coordination the entanglement between the quantum system and the reservoir near a 
phase transition can be mapped onto the entanglement which develops in an homogeneous system near a quantum 
critical point. 



II. CONCURRENCE OF THE ISING MODEL WITH GENERAL BOUNDARIES 



A. The transverse Ising model 

We start with the homogeneous, one-dimensional transverse Ising model with open boundary conditions and cou- 
pling parameter A. The two spins at the end are further connected by an additional coupling parameter k. For k = A, 
one recovers the Ising model on a ring. The full Hamiltonian thus reads 

N-l N 
i=l 1=1 

where cr^'^ are the x, z-components of the Pauli matrices. 

To solve the model one first converts all the spin matrices into spinless fermions with {ci,c|,} = Si^i' [isl [l^. 
This is done by performing a Jordan- Wigner transformation 

af = exp ^ir ^ c]c}j (c, + 4) , (2) 

af = 1 - 24c, . (3) 
An additional Bogoliubov transformation then yields (up to a constant) 

N 

n = ^ujnilirin , with (4) 



N 

Vn = ^{9n,iCi + hn,icl) , (5) 
1=1 

where gn^i, hn^i, and cj„ have to be determined numerically for arbitrary ratio k/X. Due to the unitarity of the 
Bogoliubov transformation, Eq. ([5|) is easily inverted to yield 

iv 

Cj = ^{9n,iVn + K,i-ni) ■ (6) 
n=l 

For A = 1, the energy spectrum begins at zero energy which represents the critical point. For A > 1, apart from the 
extended states at finite energies there is also an additional zero-energy "bound" state. The emergence of the bound 
state can be interpreted as a loss of coherence. Since it is connected to the appearance of a zero energy mode which 
is inherent to quantum phase transitions, we believe that this loss of coherence is a general feature that provokes the 
change in entanglement and that this view can be generalized to other systems with quantum phase transitions. For 
more details, see appendix \K\ 



B. Concurrence as information measure 



We are interested in the reduced density matrix p{i,j) represented in the basis of the eigenstates of az- It is formally 
obtained from the ground-state wave function after having integrated out all spins but the ones at position i and j. 
As measure of entanglement, we use the concurrence between the two spins, C{p{i,j)). It is defined as 



C{p{i,j)) = max{0, Ai - A2 - A3 - A4} 



(7) 
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where the are the (positive) square roots of the eigenvalues oi R = pp in descending order. The spin flipped density 
matrix is defined as p = CTj, (g) ayp*ay (8) (Jy, where the complex conjugate p* is again taken in the basis of eigenstates 
of az ■ It will be instructive to also consider the "generalized concurrence" 

C*{p{t,j)) = Xi-X2-X3-Xi . (8) 

The reduced density matrix p{i,j) p (from now on we drop the indices i and j) can be related to correlation 
functions. For this, we write the ground-state wave function as the superposition of the four states 

m = I mh^) + 1 n)\hi) + 1 imn) + 1 umi)^ (9) 

where the first ket denotes the state of the two spins at position i and j and the second ket the corresponding state 
of the rest of the spin system. The matrix element p-f|,xx = ('/'TT I'^ii)' ^-S-' ^^^^ given by Pn^u = (a^a^), where 
a± = (af±af)/2.^ 

Due to the invariance of the Hamiltonian under cr^ = ^o'f i least eight components of the reduced density matrix 
are zero (for finite N). The diagonal entries read: 

pi EE ptT^TT = (1 + i'^t) + K> + K'^|»/4 (10) 

p, EE pu^n = (1 - ('^f) + (^|> - (<^|»/4 (11) 

P3 EE pxrar = (i + ('^f) - K> - (<^|»/4 (12) 

P4 ^ Piui - (1 - «) - K> + «^|»/4 (13) 

The non-zero off-diagonal entries are 

P+ = PThll = (o-i^o-j^) (14) 

P- = PUAT = i'^t'^j) ■ (15) 

The positive square roots of the eigenvalues of R are then given by \yjp\p^ ± p+| and p2Pz ± P-|- Due to the 
semi-definiteness of the density matrix p, we can drop the absolute values, i.e., ^ p\pA ± /0+ and ^ p2Pz ± 

We now define I\ = pip^ — p2Pz = 4:{{afaj) — (o'f)(o'|)) and I2 = p\ — p^ = — (crfa-J)((Tf crj). For a homogeneous 
model, we have /i > and /2 > 0.[l^ The largest eigenvalue of Eq. ^ is thus given by Ai = ^PiP4 + 1/9+ 1 and the 
concurrence reads 

C*(z,j) = 2(|p+|-Vp5^) . (16) 

We note that the above expression also holds for the generalized boundary conditions. For a homogeneous system, it 
can be further simplified to 

r(z,j) = (0.,-l)/2 (17) 
where we introduced the total order dj = J2a=x y z \ 



C. Numerical results 



1. Open boundary conditions 

We first consider the nearest neighbor concurrence of the Ising chain with open boundaries [n = 0) for a fixed 
number of sites N = 101 as parameter of A, but for various positions relative to the end of the chain. The results 
are displayed on the left hand side of Fig. [1] As expected, the concurrence of the periodic model is approached as 
one moves inside the chain and the difference between C(50, 51) and C(i, i -f 1) of the periodic system is hardly seen. 
Nevertheless, the derivative of the concurrence with respect to the coupling parameter A, C = dC/dX, still shows 
appreciable differences for A « 1 (right hand side of Fig. [T|) . 

We also investigated the scaling behavior of the minimum of C'(l, 2), Xmin, for different systems sizes up to = 231. 
We did not find finite-size scaling behavior for the position of the minimum as is the case for the translationally 
invariant modcl[8]. The curve of C'(l,2), shown on the right hand side of Fig. [l] is thus already close to the curve 
for — > 00 with a broad minimum around Xmin ~ 1-1. 
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FIG. 1: Left hand side: The nearest neighbor concurrence of the open boundary Ising model for different locations relative to 
the end as function of A for — 101. Right hand side: The derivative of the concurrence with respect to A. 




FIG. 2: The generalized next-nearest neighbor concurrence of the open boundary Ising model for different locations relative to 
the end of the chain as function of A for = 101. 



The absence of finite-size scaling of the concurrence is also manifested in the case of the next-nearest neighbor 
concurrence for different system sizes N. Whereas for the periodic system the maximum of C{i,i + 2) decreases 
monotonically for TV oo,^ there is practically no change of C(l, 3) of the open chain for iV>51. 

In Fig. [21 the generalized next-nearest neighbor concurrence C*{i, i + 2) of the open boundary Ising model is shown 
for different locations relative to the end of the chain as function of A for N = 101. On the left hand side of Fig. [21 
results are shown for sites close to the end of the chain. Notice that the generalized concurrence becomes negative 
for i — 2,3 for A > 1 which is not related to the quantum phase transition. The crossover of the boundary behavior 
to the bulk behavior is thus discontinuous. On the right hand side of Fig. [21 the next-nearest neighbor concurrence 
approaches the result of the system with periodic boundary conditions as one moves inside the chain. 



2. Generalized boundary conditions 



We now discuss the concurrence for the generalized boundary conditions, introducing the parameter k. On the left 
hand side of Fig. [3l the generalized concurrence of the first two spins C*(l, 2) is shown as function of A for various 
coupling strengths k = 0, ..,20A and N ~ 101. For increasing k > 0, the curves indicate stronger non-analyticity at 
A ~ 1. For K>20A, the generalized concurrence becomes negative around A = 1 and is "significantly" positive only in 
the quantum limit of a strong transverse field (A < 1). A similar behavior of the concurrence is also found in the case 
of finite temperatures. [1, [l^l 
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FIG. 3: Left hand side: The generahzed nearest neighbor concurrence of the closed Ising chain for various coupling strengths 
K as function of A. Left hand side: C*(l,2). Right hand side: C(2,3). 




FIG. 4: The generalized next-nearest neighbor concurrence of the closed Ising chain for various coupling strengths k as function 
of A. Left hand side: C*(l,3). Right hand side: C*(2,4). 



On the right hand side of Fig. [31 the concurrence of the second two spins C(2, 3) is shown. All curves display similar 
behavior. There is thus a rapid crossover from the boundary to the bulk-regime and the concurrence of periodic 
boundary conditions is approached for all k as one moves further inside the chain. 

To close, we discuss the next-nearest neighbor concurrence C{i,i -\- 2) for various values of k and N = 101. On 
the left hand side of Fig. 01 the generalized concurrence of the first and the third spin, C*(l,3), is shown. For k<1, 
C*(l,3) is positive for all A. For C*(1,3) first becomes negative for A < 1. For k>1.5 C*(1,3) is negative for 

all A. On the right hand side of Fig. 21 the generalized concurrence of the second and the forth spin, C*(2,4), is 
shown. For k<X/2, the C*(2,4) is negative for A>1. For k>3A/2, the C*(2,4) is negative for A<1. Nevertheless, the 
maximum value is close to A = 1 for all cases. 

We finally note that the third neighbor concurrence remains zero for all A and all k. 



D. Summary 



To conclude, we have calculated the entanglement between qubits at the boundary of a spin chain, whose parameters 
are tuned to be near a quantum critical point. The calculations show a behavior which differs significantly from the 
that inside the bulk of the chain. Although the spins are part of the critical chain, we find no signs of the scaling 
behavior which can be found in the bulk. Still, we could identify a boundary regime, basically given by the first site, 
and a crossover regime of approximately 10 sites till the bulk behavior is reached. We use the same approach as done 
previously for the bulkj^, 1^], although it should be noted that the existence of a finite order parameter in the ordered 
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phase will change these results if the calculations were performed in the presence of an infinitesimal applied field. 

III. CONCURRENCE AT A BOUNDARY PHASE TRANSITION 

In order to observe critical behavior of the concurrence at the boundary, one has to consider a different model than 
the simple transverse Ising chain. One possibility would be to introduce an isotropic coupling from spin N to spin 1 
which would lead to an interaction term containing four fermionic operators. A simple solution is thus not possible 
anymore. In the following, we will consider a similar model, but which can easily be mapped onto the spin-boson 
model. 



A. The model 

The model with isotropic coupling between the two spins at the end is similar to the model introduced by Garst 
et al. [IJl (see also Ref. (l^l)- It describes two spin-1/2 systems attached to two different electronic reservoirs. They 
further interact among themselves through an Ising term. We can write the Hamiltonian as 

H = Hki + Hk2 + ISzlSz2 , 

k k,k' .j.L.u 

(18) 

To evaluate the concurrence of the two spins, the 4x4 reduced density matrix in the basis of the eigenstates of 5^1 
and Sz2 is needed. 

The system described by Eq. (fT8|) undergoes a Kosterlitz-Thouless transition between a phase with a doubly de- 
generate ground state and a phase with a non degenerate ground state. This transition is equivalent to that in the 
dissipative two-level systemjy, 0] as function of the strength of the dissipation. We define the dissipative two-level 
system as 

Htls = Aa, + ^ \k\blbk + Aa, ^ Vk{hl + h^) . (19) 

k k 

The strength of the dissipation can be characterized by a dimensionless parameter, a oc A^, and the model undergoes 
a transition for A — 6/lUc ^ 1, where Uc is the cutoff, and a — 1. The Kondo model can be mapped onto this model 
by taking A cx Jl and 1 — a (x Jz [Hi ■ 

To understand the equivalence between these two models, it is best to to consider the limit I/J^l (the transition 
takes place for all values of this ratio). Let us suppose that / > so that the Ising coupling is antiferromagnetic. The 
Hilbert space of the two impurities has four states. The combinations | tt) and | H) arc almost decoupled from the 
low energy states, and the transition can be analyzed by considering only the | "f |) and | 11) combinations. Thus, we 
obtain an effective two state system. The transition is driven by the spin flip processes described by the Kondo terms. 
These processes involve two simultaneous spin flips in the two reservoirs. Hence, the operator which induces these 
spin flips leads to the correspondence A ^ j\/{Iuic)- The scaling dimension of this term, in the Renormalization 
Group sense, is reduced with respect to the ordinary Kondo Hamiltonian, as two electron-hole pairs must be created. 
This implies the equivalence 2 — a Jz- Hence, the transition, which for the ordinary Kondo system takes place 
when changing the sign of Jz now requires a finite value of Jz- 

B. Calculation of the concurrence 

The 4x4 reduced density matrix can be decomposed into a 2 x 2 box involving the states | ||) and | which 
contains the matrix elements which are affected by the transition, and the remaining elements involving | ||) and 
I II) which are small, and are not modified significantly by the transition. Neglecting these couplings, we find that 
two of the four eigenvalues of the density matrix are zero. The other two are determined by the matrix 



1/1 + (a,) (a,) \ 
(a.) l-{<yz)-) 



(20) 
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where the operator a is defined using the standard notation of the dissipative two level system, Eq. 
entanglement can be written as 



IHl). The 



(21) 



The value of (ctz) is the order parameter of the transition. The value of {px), at zero temperature, can be calculated 
from 



(22) 



where £^o is the energy of the ground state. Using scaling arguments (see appendix |B|, it can be written as followed: 



A (A) _ A! 



En = < 



l-2a 

c 



UJc J 



< a < i 



< a < 1 

a - 1 
a > 1 



(23) 



where C, C" and C" are numerical constants. 

If the density matrix is calculated in the absence of a symmetry breaking field, (az) — even in the ordered phase. 
Then, from Eq. ([2T|) . the concurrence is given by C = K^x)], which is completely determined using Eqs. ([22|) and (|23|) . 
In the limit A/ujc ^ 1 the interaction with the environment strongly suppresses the entanglement. We expect unusual 
behavior of the concurrence for a = 1/2 and a = 1. The point a = 1/2 marks the loss of coherent oscillations between 
the two states [13, HI], although the ground state remains non degenerate. Following the analysis in [1], we analyze the 
behavior of dC/da, as a is the parameter which determines the position of the critical point. The strongest change 
of this quantity occurs for a = 1/2, where: 



dC_ 

da 



a=l/2 



A /UJc 
— log I -T- 

uJc V A 



(24) 



On the other hand, near a — 1 the value of dC Ida is continuous, as the influence of the critical point has a functional 
dependence, when a ttc, of the type e~'^/'-"'=~"-'. This is the standard behavior at a Kosterlitz-Thouless phase 
transition. This result suggest that the entanglement is more closely related to the presence of coherence between 
the two qubits than to the phase transition. The transition takes place well after the coherent oscillations between 
the I tl) and | |t) states are completely suppressed. We note though that with a symmetry breaking field, there is a 
discontinuity of the concurrence at the phase transition [2^ . 



IV. VON NEUMANN ENTROPY FOR DISSIPATIVE SYSTEMS 



In this section, we will use the von Neumann entropy as measure of entanglement. It is defined for any bipartite 
system with a ground-state \ip) by introducing the reduced density matrix with respect to one of the subsystem. For 
the two subsystems A and B, it reads 

E{ij) = -Tr{pA\npA) , = Tr^d V) (V-l) • (25) 

In contrast to the concurrence, an analytic expression of the reduced density matrix pA does not automatically lead 
to an analytic expression for the von Neumann entropy. In the following, we show that in the case of integrable 
dissipative models, the von Neumann entropy can be obtained analytically. We then also discuss the von Neumann 
entropy of the spin-boson model. 



A. Integrable quantum dissipative systems 

Modeling the environment by a set of harmonic oscillators 4], the canonical (integrable) model for dissipative 
systems is described by the following Hamiltonian: 
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The operators obey the canonical commutation relations which read {U = 1) 

[q,p]=i , [Xa,Pa'] = iSa,a' ■ (27) 

The coupling of the system to the bath is completely determined by the spectral function 

JffoH = 5V^<5(u;-c^„). (28) 

a 

In the following, we will consider a Ohmic bath with Jho{^) — V-^ for lu ^ lUc and Jho{^) = for lu 3> lOc, lOc being 
the cutoff frequency. 

1. Caldeira-Leggett model 

Let us first consider the free dissipative particle, i.e., we set loq 0. The model was introduced by Caldeira and 
Leggett [131 and further investigated by Hakim and Ambegaokar [5^ . The latter authors obtained the reduced density 
matrix via diagonalization of the Hamiltonian. In real space, it reads 

{x\pA\x')^e"^^^-^'^'/L , a-i^lnfl + 4) (29) 

where 77 denotes the phenomenological friction coefficient and ujc is the cutoff frequency of the bath, introduced below 
Eq. Furthermore, L ^ 00 denotes the system size and in contrast to the use of Eq. [51] in Ref. [281, here the 

normalization is crucial to assure Tip a = 1. 

In order to calculate the entropy of the system, we Taylor expand the logarithm 

n=l n=l fe=0 \ ^ 

Further we have 

MW)-^' '/fe-t(--')Vi'= (31) 

proved by induction. With the identity 
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k x/tt 



dxe^"'' (32) 



we thus obtain for the specific entropy (for general dimension d) 

S^^ [HaL^) + Men)] . (33) 

Comparing the above result with the entropy of a particle in a canonical ensemble, we identify a ^ A^^ oc T with 
A denoting the thermal de Broglie wavelength and T the temperature of the canonical ensemble. Notice that the 
entropy of a free dissipative particle shows no non-analyticity. 

2. Dissipative harmonic oscillator 

We now include the harmonic potential, i.e., ujq 7^ 0. The reduced density matrix of the damped harmonic oscillator 
is given by [Tj] 

{x\pa\x') = . &-'^(— -')^~K-+-')= ^ (34) 
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with a — and b 



g^ijy. The above expression is deduced such that the correct variance for position and 



momentum is obtained. At T = the expectation values are given by 

1 



ZUJQ TT 



UJQ 



with K ~ r\llLLiQ and 



/(«) 



1 In [(k + ^/^2— _ ^f^^~\)\ 

TT - 1 



(35) 



(36) 



The parameter k represents the friction parameter and the system experiences a crossover from coherent to incoherent 
oscillations at k = 1. 

Taylor expanding the logarithm of the entropy, Eq. (j30p . leads to the evaluation of the general n-dimensional 
integral 



da;i..(ia;„ exp I — XiAi,j'^j 



■i/2 



VdetA 



(37) 



where A is given by the translationally invariant tight-binding matrix with Ai^i — 2(a + 6), Ai^i^i = = — (a — &) 

(n + 1 = 1) and zero otherwise. The determinant of the matrix is given by its eigenvalues and reads 



detA= (2a)"(l-6/a)" ]J 



ni—1 



2b 



cos kr, 



(38) 



with km = 27rTO/n. 

The determinant can be easily evaluated for large cut-offs lUc ^ oo 29]: Considering the n-dimensional transla- 
tionally invariant, but non-Hermitian matrix Ai^i — 1, — 1 — e (n -1-1 = 1) and zero otherwise, one obtains the 
following formula: 



n 



For LOc/^o S> 1, we have 



2(1 -e) 
a/6 = 4(g2)(p2) 



— cos fc„ 



2"(1 - £)" 



n\2 



(39) 



4k . 



(l-2K^)/(K)-f — In^ 



Wo 



> 1. 



In this limit, we can thus set = Ab/a <^ 1 and the n-dimensional integral can be approximated to yield 



dx{x\p^\x) 



l-{l-ey 



with 



1 — 6^/4. Expanding the denominator as geometrical series, we then have for the entropy 

e 



S 



■ Ine + ^ ln(l - e] 



(40) 



(41) 



(42) 



In the limit e « £ ^ 1, the leading behavior of the entropy is given by 5 ~ ln(a/6). 

The determinant can also be calculated exactly as was done in Ref. [26|]. This yields the exact expression of the 
von Neumann entropy, 



S = 

2 



Ab 

a — b 



-^In 
2b 



\/a- Vb 
y/a + Vb, 



(43) 



The von Neumann entropy is plotted in Fig. [5] for various cutoff energies Uc as function of the coupling constant n. 
Notice that in all curves a crossover behavior occurs at k w 1, where coherent and incoherent oscillations interchange. 



10 




FIG. 5: The entropy S of the dissipative osciUator with Ohmic couphng as function of the dimensionless coupling strength n 
for various cut-off frequencies uuc- 



B. Spin-boson model 

In section IIIIl the spin-boson model or dissipative two-level system was already introduced and the concurrence 
was calculated in the context of a two-impurity Kondo model. Here, we want to compute the von-Neumann entropy 
for this system. Since we will discuss several bath types, the Hamiltonian without bias shall be defined with general 
coupling constants as 

^-^^.+E^^^I^fc + '^-E Y(^fc + ^I) • (44) 

fe k 

Again, the operators b^J^ resemble the bath degrees of freedom and a^, Oy, denote the Pauli spin matrices. The 
coupling constants give rise to the spectral function 

JM=E^fe'^(^-^fc)- (45) 
fe 

In the relevant low-energy regime, the spectral function is generally parameterized as a power-law, i.e., Jiyi) oc 
200;** Aq"" where a denotes the coupling constant, s the bath type (s = 1 defines the previously discussed Ohmic 
dissipation) and Ag the cutoff- frequency. The change in notation Ag) will be convenient in the context of the 

scaling approach. 

With A denoting the spin-1/2 system, the reduced density matrix of the spin-boson model is given by 

1(1 + (a.) K) \ 

Since there is no symmetry breaking field in the above Hamiltonian, we set (az) — 0. The eigenvalues are thus given 
by A± = (1 ± (crj;))/2 and the entropy reads 



2 



ln((l-K)^) /4) + (a.>ln(^i±i^ 



(47) 
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0,5 1 0,2 0,4 0.6 0.8 1 

a a 

FIG. 6: Left hand side: The entropy S of the spin-boson model with Ohmic coupling in units of ln(2) as function of the 
dimensionless couphng strength a for various cut-off frequencies Aq. The right hand side shows \n{ax) which derivative with 
respect to the couphng approaches a step-hke function in the scaling limit Aq/Aq — > 0. 



The value of (ax), at zero temperature, is given by 



where Eq is the energy of the ground-state. To obtain the ground-state energy, a scaling analysis for the free energy 
at arbitrary temperature is considered as before (see appendix [B|) . £^0(^0) and (ax) will then set the basis for our 
discussion on the entanglement properties of the spin-boson model, see Eq. (|47p . 



1. Ohmic dissipation 

In the Ohmic case js — 1), there is a phase transition at zero temperature at the critical coupling strength 
a = 1 + 0(A/Ao) (boI. l3ll|. The transition is also reflected by the renormalized tunnel matrix element Aj-en which 
reads Are„ = Ao(Ao/Ao)"/(^-") for a < 1 and A.-cn = for a > 1. 

The von Neumann entropy of the spin-boson model with Ohmic dissipation was first discussed by means of a 
renormalization group approach [s^l and later also by the thermodynamical Bethe ansatz(33|. Here, we will obtain 
the von Neumann entropy within a scaling approach which can also be extended to non-Ohmic dissipation. In this 
approach, the free energy is given by (see appendix IB|) 

With A(A) = Ao(A/Ao)", the ground state energy Eq is then given by Eq. (P5)) and the discussion is similar to the 
one in section UlI Bl For a = 1/2, we thus have 



dlTi{ax 



da 



oc In f ^ ) . (50) 
--1/2 \Ao ' 



In the scaling limit Aq/Aq — s- 00, this quantity diverges logarithmically. This is shown on the right hand side of Fig. 
[SI The entropy S of the dissipative two-level system with Ohmic coupling is plotted in Fig. [S] as function of the 
dimensionless coupling strength a for various cutoff frequencies Aq. The entropy quickly saturates after the transition 
from coherent to incoherent oscillations at a = 1/2 as can be seen in terms of \n{ax) on the left hand side of Fig. [SI 
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2. Non-Ohmic dissipation 



The calculation of _Bo(Ao) and {ax) can be extended to the spin-boson model with non-Ohmic dissipation (s ^ 1). 
In general, the dependence of the effective tunneling term on the cutoff, A(A), is: 

A(A) = Ao exp l^-i j^'' '^du)j (51) 

with the spectral function given in Eq. (|45p . A renormalized low energy term, Aj-en; can be defined by 

A,e„ = Aoe-^t''o„^'^- . (52) 

The free energy is again determined by Eq. though cannot be evaluated analytically, anymore. The scaling 

behavior of the renormalized tunneling given in Eq. (jSip is no longer a power law, as in the Ohmic case. Still, we can 
distinguish two limits: 

i) The renormalization of A(A) is slow. In this case, the integral in Eq. ((49|) is dominated by the region A ~ Aq, 
where the function in the integrand goes as A~^. The integral is dominated by its higher cutoff, Aq, and the 
contribution from the region near the lower cutoff, Aj-en, can be neglected. Then, we obtain that -F(Ao) ~ Ag/Ag. 

ii) The renormalization of A(A) is fast. In this case, the contribution to the integral in Eq. I|49p from the region 
A « Aq is small. The value of the integral is dominated by the region near A ~ Aren- As Aron is the only quantity 
with dimensions of energy needed to describe the properties of the system in this range, we expect that F(Ao) ~ Arcn- 

In the scaling limit, Aq/Aq ^ 1, the values of the two terms, Aj-en and Aq/Aq, become very different. In addition, 
there are no other energy scales which can qualitatively modify the properties of the system. We thus conclude that 
only the two terms mentioned above will contribute to the free energy. Hence, we can write: 

F(Ao) ~ max f A^en, ^\ (53) 



A_ 

In the following, we will use this conjecture to discuss super- and sub-Ohmic dissipation. 

a) Super-Ohmic dissipation. In the super-Ohmic case (s > 1), Eq. ()52|) always has a solution and, moreover, we 
can also set the lower limit of the integral to zero. This yields 

A,en = Aoe- ^o" ^"'^ « Aoe-"/(^-i^ . (54) 

For a 3> 1 we have Aicn ^ Aq, but there is no transition from localized to dclocalized behavior. 
Using Eq. (|53p in the super-Ohmic case s > 1, we can approximately write: 

(0^niax(^e-"/(^-i),^) (55) 

We thus find a transition from underdamped to overdamped oscillations at some critical coupling strength 
a~ (s-l)log(Ao/Ao). 

It is finally interesting to note that the scaling analysis discussed in Ref. is equivalent to the scheme used 
here. 

b) Sub-Ohmic dissipation. In the sub-Ohmic case (s < 1), it is not guaranteed that Eq. (|52|) has a solution. In 
general, a solution only exists when Aq/Aq is not much smaller than 1. 

The existence of a phase transition in case of a sub-Ohmic bath was first proved in Ref. [s^ . Whereas the relation 
in Eq. (j52p and a similar analysis based on flow equations for Hamiltonians j36| yields a discontinuous transition 
between the localized and delocalized regimes, detailed numerical calculations suggest that the transition is 
continuous [37| . 

Since there is a phase transition from localized to non-localized behavior, there might also be a transition 
between overdamped to underdamped oscillation. In Ref. [33|, this transition was discussed on the basis of 
spectral functions analogous to the discussion of Refs. [U, for Ohmic dissipation. It was found that for 
s > 0.5 the transition takes place for lower values of a as in the Ohmic case, e.g., for s — 0.8 and Aq/Aq = 10 the 
transition coupling strength is a* ~ 0.2. For a recent discussion on the spectral properties using the Numerical 
Renormalization Group, see Ref. ( ^40|] ) . 
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Aq/Aq 



DL(coh.) 



DL(incoh.) 



L 



a 1 



FIG. 7: Schematic picture of the different regimes in the sub-Ohmic dissipative TLS studied in the text. DL stands for the 
delocalized phase, while L denotes the localized phase. The lower blue line denotes the continuous transition studied in Ref. 
[Sj . The red line marks the boundaries of a regime characterized by a small renormalization of the tunneling rate, Eq. (I52|) . 
and coherent oscillations. 



Using Eqs. ((52|) and ([53]) yields for the sub-Olimic case the following qualitative behavior: 




delocalized regime ^ ~ 1 
localized regime ^ 1 

° An 



^" ^ (56) 



The analysis used in the previous cases leads us to expect coherent oscillations in the delocalized regime. 

We can extend the study of the sub-Ohmic case to the vicinity of the second order phase transition described in 
Ref. [ilj, which in our notation takes place for a = sAq/Aq <C 1. In this regime, which cannot be studied using 
the Franck-Condon like renormalization of Eq. (j52p , we use the renormalization scheme around the fully coherent 
state proposed in Ref. [4ll |. To one- loop order, the beta-function for the dimensionless quantity (expressed in 
our notation) k = (aA) / A then reads 

P{k) = -sk + k^. (57) 
Near the transition, in the delocalized phase, k thus scales towards zero as 



The scaling of {ax) is 



k{A) = ko[^J . (58) 



The fact that the scheme assumes a fully coherent state as a starting point implies that A is not renormalized. 
Inserting Eq. into Eq. we find: 

d{ax) , fAVAo 



OA VAo/ A2 



If we calculate {cTx) from this equation, we find that the resulting integral diverges as A — > for s < 1. This 
result implies that {ax) <C 1. For sufficiently low values of the effective cutoff, A, the value of {ax) can be 
calculated using a perturbation expansion on Aq, leading to {ax) ~ Aq/Aq. This result implies the absence of 
coherent oscillations. A schematic picture of the regimes studied for the sub-Ohmic TLS is shown in Fig. [7]. 

We finally note that the entanglement of a spin-1/2 particle coupled to a sub-Ohmic environment has recently been 
discussed in Ref. [i^ . 
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a) 



N 



b) 



Nj+1 



c) 



N, 



m-<y- 



FIG. 8: a) Measuring the system at the left end of the environment, b) Measuring part of the environment, c) Measuring the 
system and part of the environment simultaneously. 



V. NON-LOCAL INFORMATION MEASURE 



Quantum measurement is closely connected with the collapse of the wave function and due to the recent advances in 
quantum engineering, the concept of "information" has to be reconsidered when one deals with quantum mechanical 
systems. But instead of introducing a new concept of quantum information "from scratch" , one can also start with the 
measuring process and see what information can be extracted. This line was recently pursued by Zurek and coworkers 
[43I I by proposing that in a classical description, information can be obtained by measuring the environment to which 
it is coupled. This approach seems even more appropriate for quantum mechanical systems. 



A. The model and information measure 



In this section, we want to employ an information measure based on the measurement process and apply it to a 
dissipative quantum system. The model will consist of a harmonic chain with open boundaries. If the mass of the 
first (quantum mechanical) particle is large compared to the other masses, one speaks of the Rubin model [i^l, but 
for simplicity, we will choose all masses equally, here. In subsection IV Bl we will then distinguish between the spring 
constant of the bulk and of the edge. The Hamiltonian is given by 

The particle at the left end of the chain shall denote our system which is coupled to the environment B\ (filled and 
empty circles in Fig. [5^), respectively). 

a) Measuring the system means that one is only interested in the mean value of the environment. The relevant 
density matrix is thus obtained by tracing out the bath degrees of freedom. The von Neumann entropy is known 
to be a good measure to characterize the ground state. We thus have 

£;i(^) = -Tr(pilnpi) , pi = Tr^, (jV') (V-l) • (62) 

The above model can be mapped to the dissipative harmonic oscillator with Ohmic coupling. This is done by 
diagonalizing the bath modes 

— / dfc sin(fcn)x(fc) (63) 



and results in the following representation of the Hamiltonian 

2 _f / 2 2 



k ^ ^ k 
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b) In order to apply the information approach proposed in Ref. [43| , we will now pick out one of the environmental 
particles, see Fig. |5}d). Again, the dissipative system (|6ip can be brought into more familiar form by diagonalizing 
the left and the right part of the environment separately. This formally results in the problem where a quantum 
mechanical particle in a harmonic potential is coupled to two baths. But since the left and the right reservoir 
linearly couple to the same spatial coordinate, they are indistinguishable. The resulting model is thus the 
standard dissipative harmonic oscillator with modified coupling functions Xk as given by Eq. (|64|) . 

The relevant density matrix for the selected particle of the environment is now obtained by tracing out the bath 
degrees of freedom without the selected particle plus the system itself, labeled as i?2- For the von Neumann 
entropy we thus have 

£;2W --Tr(p2lnp2) , P2 = Tr^, (|^) (V-l) ■ (65) 



c) The last step is to measure both, the system at the left end of the chain and the selected particle of the bath, 
see Fig. [8};). Again, we proceed by decoupling the left and right part of the environment separately. We obtain 
the following representation of the Hamiltonian in Eq. (pT|) : 

k,i=l,2 \ / 

+ xi ^ Afc,ia;fc,i + 2^2 ^ {Kk,iXk,i + Kk,2Xk,2) (66) 

fc k 

In contrary to case b) , here there is a distinction between the two resulting non- interacting reservoirs since one 
bath is coupled to two particles whereas the other bath only couples to the environmental particle. There is no 
way of preforming a unitary transformation such that the two reservoirs act as one. 

A similar type of problem has been analyzed by Kohler and Sols where two different baths were coupled to the 
momentum and to the spatial coordinate, respectively Also from the two-channel Kondo model it is known 
that two baths can significantly alter the system behavior due to the simultaneous measurement process. [4^ 
We thus expect that effects of quantum frustration are contained in the employed information measure. 

The von Neumann entropy of the subsystem is given by tracing out the degrees of freedom of the remaining 
bath B3 

E3W = -Tripsin ps) , ps^TreMm • (67) 

The measure of information which is contained by measuring parts of the environment as proposed by Ref. [43] is 
now given by 

Eii;) ^Es{^)-E^ (V-) - E2 (^) . (68) 

In the following, we will set N2 = 0, i.e., we investigate the entanglement between the two ends. In the context of 
spin-models, the long-distance entanglement was recently considered using as measure of entanglement the concurrence 
[47| . On the other hand, it was shown that the above measure based on the von Neumann entropy only captures 
classical correlations if it is positive [il] . 



B. Entanglement between the two ends 

For explicit calculations, we will consider a simplified version of the above model and neglect the reservoir to the 
right, i.e., we will set N2 ^ Q va Fig. [S]b), c). This amounts to the following question: What is the entanglement 
between the two ends of a harmonic chain as function of the system size N . 

The chain is confined by the masses at xq and xn- The diagonalization of the harmonic chain for finite length 
N -I yields 

N-l 

Xn , ^ y ^ S\n[krnn)Xm , Ki = • (69) 

^ rn— 1 

The eigenvalues are given by = 4/^ sin^(fcm/2). Here, we have introduced an extra spring constant fs for the 
masses of the bath to contrast it from the spring constant that connects the two masses at the end with the chain, 
denoted by /. 
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In the following, we will neglect finite size effects and only consider the case where there are two particles at the end. 
The case of one particle is then simply obtained by neglecting the second particle and the transformed Hamiltonian 
reads 

^ - E y + ^'^^^ + E (f + - E -™-™'^^) ' (70) 

i=l,2 m 1=1,2 



zl = ^ sin(A:„) , cl = (-1)"^ MK.) (71) 



with {N N) 



and gj=i = xq and qi=2 ^ xn- 

To obtain the von Neumann entropy of the various subsystems listed in a)-c) , we first need to compute the reduced 
density matrix. The reduced density matrix of dissipative systems is commonly represented as a path integral where 
the bath degrees of freedom have been integrated out [1, 0| : 



{q"\pA\q') (X / VqeM'SPiq] - ^^Jg]) (72) 



Here Sg denotes the Euclidean action of the system and S^^ji the influence on the system due to the environment. 
For one particle coupled on a linear chain with coupling constant Cm, we have 



with the Fourier transform 



9(t) = ^ ^ g„ exp(w„r) , i/„ = 2?™//? . (74) 



For two particles coupled to both ends of a linear chain with coupling constants Cm^'^\ we have 

oiE) r 1 kmgn + Cmg»P fr,f■^ 



n ' m 



Notice that there is no potential renormalization in our model originating from the harmonic chain. 

Since we have already discussed the von Neumann entropy for a dissipative particle in a harmonic potential, we 
are left with the case of two particle, see Fig. c). With the coupling coefficient of Eq. ([7T|) . the effective action 

S — Sg + S\^ji can be written as S* = 5*1 + 5*2 + Sint with 

E \\T.^^l+^l-^li^n))Wn? (76) 
i=l,2 ^ n 

the effective action of particle i coupled to the dissipative environment and 

S^nt = ^]^Y.^]{v-^)^e{qlql) {77) 

n 

the effective action describing the interaction between the two particles through the environment. In the above 
equations, we further defined the potential renormalization LOr{i^n) and the (system-size independent) effective splitting 
parameter ojiivn) as 

, ,21,, \ 1 /I siv?{km) , . 

^^"^"^„t'r-^ + 4/.sin^(W2) ' ^''^ 

u;}{vn) = hir — — . (79) 
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By a unitary transformation, = {qj^ ± qf-^)/\/2, the two modes can be decoupled, i.e., 5' = S-^- + 

^^E(^« + ^o-^'K)±^?K)/^)knP (80) 

1=1,2 ^ n 

The physical behavior of dissipative models is determined by the low-frequency modes of the bath. The action 
can thus be interpreted as the action of two harmonic oscillators with the effective frequencies uj± — -^tD^ ± (coi/N)^ 
where lo = \/ujq — uj^- For i/„ — > 0, we further have cu'^ = ujj = /b/2. 

For the chain with equal spring constant /b = f, we have ujq = //2 = /b/2 and thus a) = 0, which indicates a 
phase-transition to a localized state. For f > fs, we can use the results of the entropy of an harmonic oscillator. In 
the expression of the entropy Eq. (|43|) . only the combination {q^){p'^) enters, such that the only dependence on the 
system size N is contained in the term 

(.^>(p^>->^ln^ . (81) 

TT UJ± 

Expanding the logarithm, the linear term cancels and we thus have for the information measure for two particles at 
the end of a harmonic chain with length N the following scaling behavior: 

I = E{A + B)-E{A)-E{B)^^ (82) 



VI. SUMMARY 



In this article, we have investigated the entanglement of quantum systems at the boundary. We have first calculated 
the entanglement between qubits at the boundary of a spin chain, whose parameters are tuned to be near a quantum 
critical point. The calculations show a behavior which significantly differs from that inside the bulk of the chain. 
Although the spins are part of the critical chain, we find no signs of the scaling behavior which can be found in the 
bulk. We use the same approach as done previously for bulk spins although it should be noted that the existence 
of a finite order parameter in the ordered phase will change these results if the calculations are performed in the 
presence of an infinitesimal applied field. 

We have also considered the entanglement between two Ising-coupled spins connected to a dissipative environment 
and which undergo a local quantum phase transition. The system which we have studied belongs to the generic class 
of systems with a Kosterlitz-Thouless transition at zero temperature, like the Kondo model or the dissipative two level 
system. The most remarkable feature of our results is that the entanglement properties show a pronounced change 
at the parameter values where the coherent quantum oscillations between the qubits are lost. 

In the second part of this article, the entanglement properties of dissipative systems were investigated using the von 
Neumann entropy. We first discuss two integrable dissipative quantum systems - the free dissipative particle and the 
dissipative harmonic oscillator - and calculated the von Neumann entropy. In the former case, we found an analogy 
to the entropy of a canonical ensemble at temperature T . The case of the harmonic oscillator is the more interesting 
one since it exhibits a transition from underdamped to overdamped oscillations. This transition is also manifested 
in the entropy, but not as strongly as in the case of the spin-boson model. This is probably due to the absence of a 
quantum phase transition and that the model can be adequately treated by semi-classical methods as done e.g. in 
the context of the fluctuation-dissipation theorem . 

We also calculated the von Neumann entropy for the spin-boson model on the basis of a scaling approach for the 
free energy. Only in the Ohmic case, the resulting integral, i.e., the ground-state energy, could be evaluated and we 
analyzed the behavior at the transition from underdamped to overdamped oscillations. We found that the change of 
the logarithm of {ax) with respect to the coupling strength a is strongly pronounced at the Toulouse point. In the 
non-Ohmic case, we argued that the crossover between coherent and decoherent oscillation takes place when the value 
of ((Tj;) becomes comparable to the result obtained using a perturbation expansion in the tunneling matrix (as it is 
the case for Ohmic dissipation). In this framework, we can also discuss the super-Ohmic and sub-Ohmic dissipative 
two-level system, respectively. We conclude that entanglement properties are closely connected to the transition of 
coherent to incoherent tunneling. 

In the third part of this paper, we have applied an extended measure of quantum information to a simple model, 
describing a chain of harmonically coupled particles. We argued that this measure can be applied to relate particles 
of arbitrary distance (or arbitrary regions of the chain) and that it incorporates features of quantum frustration. We 
calculated explicitly the information measure which relates the two particles at the two ends of the harmonic chain 
which decays algebraically with the system size. 



18 



VII. ACKNOWLEDGMENTS 

Funding from FCT (Portugal) grant PTDC/FIS/64404/2006 and from MFC (Spain) grant FIS2004-06490-C03-01 
is acknowledged. 

APPENDIX A: JORDAN- WIGNER AND BOGOLJUBOV TRANSFORMATION 

In this appendix, we start the discussion with the slightly more general anisotropic spin-1/2 Heisenberg model in a 
homogeneous magnetic field, which is given by 



(Al) 



where sf /2, cr" denoting the Pauli matrices with a = x,y,z. 

Introducing the new operators Oi = sf + is^ (which leads to sf = a|ai — 1/2), one now performs a Jordan- Wigner 
transformation 1491 



Ci = exp < in OjUj > , a,; = exp < —in CjCj > Ci 



(A2) 



With now anti-commuting c^-operators {ci, c^j} — the Hamiltonian can thus be written in Fermion operators as 
H = E + /».c.)^^^(cl4i + h.c.) + - l/2)(4,c,+i - 1/2) - HcU - 1/2)1 • (A3) 

For the Ising model in a transverse field we set Jy ^ Jz ^ ^ which yields 

i i 

where we chose fixed boundary conditions since aj^ai ^ c^jqCi. For general boundary conditions, i.e., k 7^ in Eq. ([1]), 
we will neglect the boundary term that involves the operator exp(i7r A^i) i'^ order to preserve the bilinearity of 
the model, see Ref. ITsl. Using a more general notation for a bilinear Hamiltonian 



(A5) 



the formal diagonalizing of the above Hamiltonian via a Bogoljubov transformation '50'| leads to 



N 1 



N 



N 



(A6) 



where rja = ^^gaACi + ha.icl are the operators in the diagonal basis. 

The Bogoljubov transformation, i.e., the determination of the new eigenenergies as well as the new operators 
through ga^i and ha,i, is equivalent to solving the eigenvalue problem of the matrix {A — B){A + B). For the Ising 
model with open boundary condition, this is equivalent to the problem of a one-dimcnsional chain with an impurity 
at the first site, i.e., {A - B){A + B) ^ {Jx/'^fHrB with 



Af-l 



JV 



Htb = - 51 MAu+i + h.c.) + ^i''t\tl + (/^2 + 4) ^ t\u 



(A7) 



i=2 



where /i = Ah/ Jx. 

We now want to analyze the eigenvectors and eigenenergies of the tight-binding model. For 1 < n < A^, the 
eigenvectors Xi and eigenenergies = tt'^/( Ji,/4)^ are given by the following equations: 



-2nXn-l + {^l^ + A- Cb'l)Xn - 2^Xn+l = 



(A8) 
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With the Ansatz a;„ 



we have for the extended states 



4 — cos ko 



(A9) 



We are interested in the hmit N oo where boundary conditions can be disregarded since k will be continuous. For 
> 0, we then have cumax = 2 + ^, Cormn = |2 - m|. 

For ji = 2, the continuum starts at zero energy which represents the critical point. For fi < 2, there is an additional 
"bound" state, i.e., x„ ^ e' 



and for n = 1, we have J)^ = fi 



. For n > 1, we then have 

cl^g = ^2 ^ 4 - 2^(6" + e-'^) , 
^ — 2fie^'^. This leads to the solution e^" = fi/2 and thus (for N — > oo) 

^ . 



(AlO) 



(All) 



The restriction fi < 2 follows from the condition k > 0, i.e., a normalizable eigenfunction. The emergence of the 
bound state can be interpreted as a loss of coherence. Since it is connected to the appearance of a zero energy mode 
which is inherent to a quantum phase transition, we believe that this view point can be generalized to other quantum 
phase transitions. 

With i J2i ^ii — —hN/2, we have for the ground-state energy 



(A12) 



With the HcUmann-Feynman theorem (m^) = —j^dhEo, we have for N ~* oo 

h — ^Jx cos k 
' J2/4 + cos fc 



{m^) = I dk 



(A13) 



At the critical point h = he = J/2, this leads to the logarithmic divergence of dh{inz 
For finite temperatures, we have 



F : 



~kT 



^kTN I dk\n[2cosh{(3uj{k)/2)] 

10 



(AM) 



With (my 



jfdhF, we have 



dk 



Jx COS k 



7o Jj/4 + /i2 -/ijj. COS fc 



tanh 



\/j^/4 + /i^ - hJx cos k 



(A15) 



The singularity of dh{mz) at h = he = J/2 is thus suppressed for T > 0. 



APPENDIX B: CALCULATION OF THE FREE ENERGY OF THE DISSIPATIVE TLS 



We calculate the free ener gy o f the dissipative two level system following the scaling approach discussed for the 
Kondo problem in Refs. (Sll. [53|. and formulated in a more general way in Ref. [ssj . For the general long-ranged 
Ising model, the scaling approach was first applied by Kosterlitz [34 1. 

The partition function of the model can be expanded in powers of as 



^ '•'0 ■'0 A A 10„ 



(Bl) 



ij — l,..,2n 



where /[(t^ — "JjO/tc] denotes the interaction between the kinks located at positions Ti and Tj. A term in the series 
is schematically depicted in Fig. [5]. The scaling procedure lowers the short time cutoff of the theory from Tc to 
Tc — dTc- This process removes from each term in the sum in Eq. (jBl[) details at times shorter than Tc — dTc- The 
rescaling Tc — > Tc — dr^ imphes the change A A(l -I- dTc/Tc). The dependence of /[(r, — Tj)/Tc] leads to another 
rescaling, which can be included in a global renormalization of A [sil WX [ssj . In addition, configurations with an 
instanton-antiinstanton pair at distances between Tc and Tc — dTc have to be replaced by configurations where this 
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Q Xi ^2 % 'C* (3 

FIG. 9: Sketch of the instanton pairs which renormalizes the calculation of the free energy of the dissipative TLS. 



pair is absent, as schematically shown in Fig. [S]. The number of removed pairs is proportional to cLtc/tc. The center 
of the pair can be anywhere in the interval < r < /3. The final effect is the rescaling: 



Z ^ Z{1 + A'^f3dT^) 

Writing Z as Z — e~^^ , where F is the free energy, Eq. (|B2p can be written as: 

dF 



In the Ohmic case, the dependence of A on Tc = is 

A(A) = Ao 

and, finally, we find the following relation: 



A^(re 



dF 


rA(A)i 


2 








A 









2q-2 



This equation ceases to be valid for A ~ Aien- For finite temperatures, we obtain 

Mo gp 
It M' 



F(T) = / —dA. 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



It is interesting to apply this analysis to a free two level system. The value of Aq does not change under scaling. 
We find the following expression: 



OA 1 A < Ao 



Inserting this expression into Eq. (|B6p . we obtain 



F(T) 



^ Ao<r 

Ao T«Ao 



and, finally: 

^ ^ dAo ^ \ l T < Ao 
in qualitative agreement with the exact result (a^) = tanh(Ao/T) 



dF r Ao « T 
a^^ — = < ^ 



(B7) 



(B8) 



(B9) 
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